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We investigate to what extent the oscillation or conversion of neutrinos en- 
hances the expected event rate of the supernova relic neutrino background (SRN) 
at the Super-Kamiokande detector (SK). The SRN i? e 's can be almost completely 
| exchanged with i^-like neutrinos by the MSW oscillation under the inverse mass 

hierarchy with Am 2 ~ 10 _8 -10 5 [eV 2 ], or by the magnetic moment of Majorana 
neutrinos with [i v > 10~ 12 /xs and Am 2 ~ 10 -4 -10° [eV 2 ]. In the standard calcu- 
lation of the SRN flux, the event rate of the SRN at the SK in the observable 
energy range of 15-40 MeV can be enhanced from 1. 2 yr" 1 to 2.4 yr _1 if all v^s are 
exchanged with i/^-like neutrinos. The enhancement is prominent especially in the 
high energy range ( > 25 MeV). In the astrophysically optimistic calculation, the 
event rate becomes as high as 9.4 yr -1 . Because the theoretical upper bound of the 
\^ , SRN events without oscillation is about 5 yr -1 taking account of the various astro- 

0^ ' physical uncertainties, we might have to resort to the neutrino oscillation if more 

than 5 events in a year, as well as a significantly harder spectrum, were observed in 
the SK. 
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There are two astrophysical neutrino sources which have already been detected until present: 
the Sun and the Supernova 1987A, and the source which is expected to become the third by the 
observation of the Super-Kamiokande detector (SK) [EJ is the supernova relic neutrino background 
(SRN). The SRN is the accumulation of the neutrinos emitted from past supernovae which have ever 



occurred in the universe |2|-|l 1| . The detection of the SRN would provide us valuable information 
on the history of the supernova rate and hence the evolution of galaxies, as well as the properties 
of the supernova neutrinos. In the earlier paper pj{ |, we calculated the SRN flux and spectrum by 
using the standard model of the galaxy evolution and predicted the realistic event rate of the SRN 
u e 's at the SK. The expected event rate is unfortunately very low: 1.2 yr -1 as the central value 
and 4.7 yr -1 as the optimistic value considering the astrophysical uncertainties; both are the rate 
in the observable energy range of recoil positrons, 15-40 MeV. The supernova rate is much higher 
in the early phase of the galaxy evolution than in the present, but neutrinos emitted from such old 
supernovae are considerably redshifted in energy and fall below the observable energy range. 

However, neutrino oscillation or conversion of the supernova neutrinos may enhance the expected 
event rate, and such effects were not considered in the previous paper [^,[0]]. We consider here 
to what extent the SRN events can be enhanced by the neutrino oscillation or conversion. The 
observable energy window is relatively high compared with the average energy of supernova neutri- 
nos, and the reaction cross section of the detection (p e p — > e + n) is proportional to E^, where E v is 
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the neutrino energy. Therefore, if P e 's are exchanged with z/^P^'s or z/ T (p T )'s, the SRN event rate 
becomes larger because of the higher temperature (or average energy) of f M -like neutrinos than P e 's 
|l2f . (In collapse-driven supernovae, z^'s, v T 's and their antiparticles can approximately be treated 
as identical particles.) The first of such possibilities is the vacuum oscillation of P e <-> P p , which 
occurs if the vacuum mixing angle is quite large compared with that of the quark sector. Although 
such large mixing is suggested from some of the candidates for the solution of the solar neutrino 
problem, such as the 'just-so' solution Jl3| or the large-angle MSW solution |Q, the conversion 
probability of P e 's and P^'s would at most be 1/2, with the maximum vacuum mixing, after the 
phase averaging. On the other hand, more efficient conversion process exists, that is, the resonant 
oscillation or conversion induced by the effective matter potential in a star (well known as the MSW 
effect |15| ). However, the ordinary MSW oscillation induced by the mixing of the flavor eigenstates 
and the mass eigenstates is relevant only to neutrinos and not to antineutrinos, under the direct 
mass hierarchy of neutrinos. The most easily detectable flavor, P e , is not exchanged with P M by this 
ordinary MSW effect. However, there exist still some possibilities of almost complete conversion 
of P e 's and z/^-like neutrinos: one is the MSW effect under the inverse mass hierarchy, and another 
is the resonant spin-flavor conversion of P e «-> za, induced by the flavor changing magnetic moment 
of Major ana neutrinos. 

There are actually some particle-physics models which give the inverse mass hierarchy jl(i| |l8[ , 
and some possible phenomenological consequences have also been discussed Jl^j2^]. The MSW 
resonance condition, 

^f C0S 29=V2G F Y e ^- , (1) 

can be satisfied with Am 2 < 10 5 [eV 2 ] above the neutrino sphere in supernovae, where Am 2 = 
m?_ — m% (defined as positive under the inverse mass hierarchy), 9 is the vacuum mixing angle, 
Y e the electron number fraction per nucleon, m p the proton mass, and p the density. We use 
here p ~ 10 11 [g/cm 3 ] and Y e ~ 0.37 at the neutrino sphere pi]], and also use E v ~ 20 MeV and 
cos 29 ~ 1. The another condition necessary for the conversion, i.e., the adiabaticity condtion is 
given by 
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at the resonance layer, where r is the radius from the center of the star, and n e the number density 
of electrons. This condition is satisfied around the neutrino sphere when sin 2 29 > 3 x 10~ 9 , where 
we assume the right-hand side of Eq. (|J) is ~ (30 km) -1 , i.e., the typical radius of the neutrino 
sphere. When Am 2 is smaller than ~ 10 5 [eV 2 ], the resonance occurs in more outer region in 
the star, and necessary mixing, sin 2 29 for the adiabatic condition becomes larger. In order to 
show this, we give here a rough estimate of the adiabaticity. Let us assume the density profile 
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is described by a power-law: p oc r . If we neglect the change in Y e , the resonance layer, 
scales as r~^ s oc Am 2 , from Eq. (|l|). On the other hand, the right-hand side of the adiabaticity 
condition, Eq. (||), scales as oc r -1 . Therefore, from the relation of r res and Am 2 , the right-hand 
side of Eq. (||) scales as oc (Am 2 ) 11 ". Then it is easy to see that the necessary mixing, sin 2 29 
scales as oc (Am 2 ) 1 /* 1-1 , i.e., 
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The numerical calculation |^2| shows that the value of k lies in the range of 2-3, and hence the 
lower bound for Am 2 to keep the right-hand side of Eq. (||) below ~ 1 is Am 2 > 2 x 10~ 8 [eV 2 ] 
(fc = 3) or 10~ 12 [eV 2 ] (k = 2). Therefore the mass range relevant to the MSW oscillation in 
supernovae is Am 2 ~ 10 -8 -10 5 [eV 2 ], and with these values of the parameters, the conversion of 
VeS and u^'s in the SRN spectrum is possible. 

Also for the magnetic moment of neutrinos, there are some practical particle-physics models 
fj3| , p4f which give the magnetic moment of about \l v ~ 10 -12 {1b, and the phenomenological con- 
sequences of the resonant conversion in the Sun or supernovae have been discussed p5|-|3"l|]. The 
conversion of D e 's and i/„'s (or i/ r 's) is especially effective in the isotopically neutral region of mas- 
sive stars, i.e., above the iron core and below the hydrogen envelope For supernovae with the 
solar metallicity, the relevant mass range jn) is Am 2 ~ 10 _4 -10 [eV 2 ], and necessary magnetic 
interaction is /ip e j, ^ 10 _12 (10 9 G/Bo)^s, where B is the magnetic field at the surface of the iron 
core and 10 9 [Gauss] is an reasonable value inferred from the observed magnetic fields of white 
dwarfs fs^] . (Conversion may occur also with Am 2 larger than 10° [cV 2 ], but this corresponds 
to the resonance in the collapsed iron core, and the detailed calculation is quite difficult.) The 
metallicity of the progenitor affects the conversion probability |3l[] , and the conversion can occur 
also with smaller Am 2 than the above range in lower metallicity stars. Such effects may be impor- 
tant in supernovae in the early phase of galaxy evolution. However, it should be noted that the 
magnetic moment of neutrinos interacts only with the transverse magnetic fields. Considering the 
global structure of magnetic fields in the collapsing star, the conversion probability changes with 
the direction of neutrinos. The conversion probability after averaging over direction, however, can 
be close to the unity, with some appropriate parameters. In Fig. [IJ the examples of the evolution 
of the conversion probability due to this mechanism are shown, by using the precollapse model of 
Nomoto and Hashimoto (1988) ||. 

Considering these possibilities, we calculate the SRN flux and event rate at the SK under the 
condition of the complete conversion, i.e., P = 1, where P is the conversion probability of V e 's 
and i/ M -like neutrinos. Because the SRN is the linear superposition of the (redshifted) spectrum of 
supernova neutrinos, it is easy to get the SRN spectrum with any values of P: 

F{E V ) = (1 - P)F Ve {E v ) + PF Vii {E v ) , (4) 

where F is the differential flux of the SRN and E v is the neutrino energy. We assume that P 
is independent of E v throughout this paper. In fact, in the both mechanisms described above, 
the conversion probability generally depends on E v . Because only the neutrinos with high energy 
(> 15 MeV) can be observed by the SK, the conversion in the whole energy range is not necessary 
to enhance the event rate of the SK, but the conversion above 15 MeV is sufficient. 

It should also be noted that strong conversion of v^s and z/„-like neutrinos is disfavored by the 
data of neutrinos from SN1987A, observed in Kamiokande II and 1MB J35|. Smirnov, Spergel, 
and Bahcall |56| concluded that P > 0.35 is excluded by the 99 % confidence level. Although 
Kernan and Krauss |3Tj arrived at the opposite conclusion that a significant oscillation is favored 
by the data, the obtained parameters of supernova neutrino spectrum are quite different from 
the theoretically plausible parameters. Jegerlehner, Neubig, and Raffelt J3j| also concluded by 



3 



using the maximum likelihood method that the fitted temperature of SN1987A data should be 
less than 4.7 MeV with the confidence level of 95.4 %, which is considerably small compared with 
the theoretically predicted temperatures. However, because of the statistical uncertainties, 
one cannot completely exclude the possibility of the full conversion. Actually, the above values 
of the confidence level are determined by the Kolmogolov-Smirnov test in Ref. |3(| and by the 
difference of the value of the likelihood function in Ref. pg[ ], both of the two methods suffer 
considerable uncertainty when the sample size is small. Therefore, we consider the case of the 
complete conversion of P e 's and i/^-like neutrinos, as the simplest model. 

In the earlier calculation which did not take account of the possibility of neutrino oscillation 
|TT[ , the spectrum of the neutrino emission from a supernova is assumed to obey the Fermi-Dirac 
distribution with zero chemical potential |12|. All stars in the mass range of 8-60 M Q are assumed 
to end their life with a collapse-driven supernova, and the region is divided into the three sub- 
regions: 8-12.5, 12.5-20, and 20-60 Mq. The total energy {Ev e ) and temperature (Tp e ) of i> e 's 
emitted from a supernova in these three mass regions are represented by those of the supernovae 
with the mass of 10, 15, and 25 Mq, respectively, by using the results of Woosley, Wilson, and 
Mayle (7). Their calculation gives (E Pe , T P J = (4.8, 4.0), (6.0, 5.0), and (11, 5.3) in units of 
(10 52 ergs, MeV), for 10, 15, and 25 M Q stars, respectively. In addition to this modeling, here 
we assume that the total energy of neutrinos is the same for v e 's and f M -like neutrinos while the 
temperature of v^-like neutrinos is higher than T Pe by a factor of 1.4 for all supernovae, in order 
to include the effect of neutrino oscillation. This simple model is inferred from the simulation of 
neutrino emission from a supernova of about 20 M Q up to 18 [sec] after the collapse by Wilson 
and Mayle [[39J (see also Ref. @|). Following their calculation, the total energy and the average 
energy of each neutrino are (4.6, 15.3) for i? e 's and (4.9, 21.7) for v^-like neutrinos, again in units 
of (10 52 ergs, MeV). With these assumptions, we can calculate the SRN flux and spectrum in the 
case that all or a part of P e 's are exchanged with ^-like neutrinos. Fig. |2| shows the calculated 
SRN differential number flux. The thick-solid line is the standard SRN spectrum obtained in the 
previous work Oj] without any neutrino oscillation or conversion. The thick-dashed line is also 
the standard SRN spectrum, but when all P e 's are completely exchanged with i^-like neutrinos 
(P = 1). The thick-dotted lines are the same with the thick-dashed line, but with different values 
of T v jTy^ ratio: 1.2 (lower line) and 1.6 (upper line). By these lines, one can estimate how the 
uncertainty of T v /Tp e ratio affects the SRN spectrum. The thin-solid line is the 'optimistic' SRN 
spectrum (no conversion), with a larger value of the Hubble constant and luminosity density of 
galaxies, and also using a different model of galaxy evolution which enhances the SRN flux in high 
energy region (see Ref. jll|). The thin-dashed line is the same with the thin-solid line, but with 
the condition of P = 1. As expected, the neutrino oscillation increases the neutrino energy and 
the spectrum becomes harder. This hardening is expected to enhance the event rate at the SK 
detector. 

Fig. H shows the expected event rate at the SK as a function of kinetic energy of recoil positrons, 
which are produced by the reaction v e p — > ne + . All the lines correspond to those of the SRN flux 
shown in Fig. p| We use 9.72 xl0~ i4 E e p e cm 2 as the cross section of the v e absorption reaction, 
where E e and p e are the energy and momentum of recoil positrons, measured in MeV. The fiducial 
volume of the SK for the SRN observation is 22,000 tons (same with the solar neutrino observation), 
and the detection efficiency is also taken into account, which is provided by Totsuka jy]. It is 
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apparent from this figure that the event rate is enhanced by the neutrino oscillation, especially in 
the high energy region ( > 25 MeV). In Table ||, the energy-integrated event rate in the observable 
energy range of 15-40 MeV is shown for the models calculated in this paper. With T v /Ta e = 1.4, 
the event rate of the standard calculation is enhanced up to 2.4 yr~* from 1.2, and the optimistic 
calculation up to 9.4 yr _1 from 4.7, if the conversion is complete. The important property of the 
effect of oscillation is that the enhancement of the event rate is prominent in the higher energy 
range. Therefore, the event rate in the range of 25-40 MeV is also shown in the table. About four 
events in a year are possible in this high energy range if we use the optimistic model. 

In conclusion, the optimistic event rate of 4.7 yr _1 can be raised up to 9.4 yr by the complete 
neutrino oscillation or conversion of rVs and f M -like neutrinos. Because the rate 4.7 yr -1 in 15- 
40 MeV is considered to be a theoretical upper bound taking account of various astrophysical 
uncertainties, we might have to resort to the neutrino oscillation or conversion if more than 5 
events were observed in the SK detector in a year. In this case, the events in the high energy 
range ( > 25 MeV) will be especially enhanced compared with the event rate spectrum without 
any oscillation or conversion. The theoretical upper bound on the SRN events in the SK becomes 
9.4 yr -1 when the neutrino oscillation is taken into account, but since the complete conversion 
seems rather unlikely considering the SN1987A data, this upper bound should be considered to be 
quite conservative one. 
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TABLE I. Expected Event Rate at the Super-Kamiokande Detector 







Event Rate [yr x ] 




Astrophysical Model 




15-40 MeV 


25-40 MeV 


Standard 


1 


1.2 


0.35 




1.2 


1.8 


0.66 




1.4 


2.4 


1.0 




1.6 


2.9 


1.35 


Optimistic 


1 


4.7 


1.42 




1.4 


9.4 


4.0 



Note. — The conversion probability, P, is assumed to be the unity, i.e., the complete conversion of 
D e 's and i/^-like neutrinos. The energy ranges are those of the kinetic energy of recoil positrons. 
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FIG. 1. The evolution of the conversion probability of v e <-> in the isotopically neutral region of 
a massive star, due to the magnetic moment of Majorana neutrinos. The 4 Mq helium core model of 
Nomoto and Hashimoto is used, and ji v = 10~ 12 ^b is assumed. The strength of magnetic fields is 
assumed to be Bo = 5 x 10 9 [Gauss], where Bo is B at the surface of the iron core, and a global magnetic 
dipole is also assumed. The values of Am 2 / E v are 3 x 10" 4 (solid line), 1 x 10~ 4 (dashed line), 5 x 1(T 5 
(dot-short-dashed line), and 1 x 10 -5 (dot-long-dashed line), in units of [eV 2 /MeV]. 
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FIG. 2. The differential number flux of the SRN as a function of neutrino energy. The thick-solid line 
shows the standard spectrum without any oscillation or conversion of neutrinos. The thick-dashed line is 
the SRN spectrum when the original zVs are completely exchanged with i^-like neutrinos. The ratio of 
temperature, T v /T& e is assumed to be 1.4, and the thick-dotted lines are the spectra with T Vfi /Tp s — 1.6 
(upper) and 1.2 (lower). The thin-solid line is the optimistic SRN flux considering the astrophysical 
uncertainties (no conversion of neutrinos). The thin-dashed line is the SRN flux using the optimistic 
model, and with the complete conversion of v e 's and v^-like neutrinos, again assuming Tv^/Ti?^ = 1.4. 
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FIG. 3. The expected event rate of the SRN v e 's at the Super-Kamiokande detector, as a function of 
the kinetic energy of recoil positrons. All the lines correspond to those in Fig. ^, in which the SRN flux 
is shown. Note that the observable energy range is about 15-40 MeV because of the other background 
events. 
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